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Abstract. Non relativistic quark models have been invoked to support the state- 
ment of a 1/rriQ duality violation in semileptonic B decays. However, as we recall, we 
have very explicitly shown that no 1/ tuq duality violation is present in totally integrated 
rates, in such quark models. Then : 1) it is shown that some contributions under dis- 
cussion are misleadingly described as "1/mg duality violation" ; as they stand, they are 
in fact parametrically much smaller : they are (^(l/mg) 5 / 2 ) when properly referred to 
the total rate under discussion ; numerically this is below 10~ 2 ; 2) moreover, relying on 
our rigourous expansion of the harmonic oscillator model around the Shifman-Voloshin 
limit, it is shown that even such small terms are not present at all in the total rate, and 
must therefore merely cancel ; 3) finally, for physical masses, by an exact treatment of 
this particular, but not too unrealistic model, one finds a departure from free-quark decay 
rate only at the 10~ 3 level ; moreover, this departure is totally explained by the OPE 
to a high precision (with an error of 10~ 5 only), provided we include the l/(mg) 3 terms. 
Further terms are strongly suppressed. On the whole, from the various detailed studies 
which have been made, it appears highly probable that equality between inclusive decay 
and OPE holds very well for total decay rates in QCD, even for a realistic spectrum and 
kinematics, provided the OPE expansion itself is pushed sufficiently far. In other words, 
not only 1/mg duality violation cannot be present, but possible duality violation, at high 
orders, is small. 
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1 Setting the controversy 

The Minnesota group [2] has established and discussed in detail the question of duality in 
QCD [2], as well as in the QCD 2 model [3] . In particular, they have extensively argued 
the famous result that there is no 1/mg term in the ratio of inclusive to free quark total 
decay rate, in semi-leptonic decays of heavy quarks, a result that was already presented in 
ref. [4]. This statement derives from the OPE. "Duality" means equality of inclusive 
rate with OPE, and, in this OPE, no 1/mg term is present. Note therefore that, to this 
order included, the inclusive rate should equate the free quark one. On the contrary, in 
ref. [5], Nathan Isgur has pointed out effects presented as a "1/mg duality violation" , a 
statement which has been often referred to. This would hold in quark models and extend 
presumably to QCD, since it has a very general structure. 

In the papers [9, 10, 11], we have studied the question for totally inclusive rates, 
using the same type of hadronic models as used by Isgur, namely quark models ; more 
precisely, we use rigourously treated, strictly non relativistic quark models, which have the 
advantage that one can calculate directly the inclusive rates as sums over bound states. 
Then, we do not find duality violation at this order 1/mg in the integrated rates, 
i.e. the rate is found to be in accordance with the OPE expectation. An important fact 
related to this conclusion is that the Bjorken and Voloshin sum rules are automatically 
satisfied in this model. We have explained here where our disagreement with ref. [5] is 
precisely taking place in ref. [1] . 

The possibility suggested in ref. [5] is continuing to be invoked in several experi- 
mental or theoretical papers in spite of the many proofs which have been provided to 
the contrary, and although no indication has been found of such an effect, as admitted 
in these papers themselves. Then, we think it worth returning to the question in de- 
tail. This paper extends our contribution to the CERN Yellow book [1]. The present 
paper contains much more complete discussions. In particular, additional informations 
on complementary aspects often intervening in the discussions are included (sections 4 
and 7). 

2 1/rriQ expansions in the non relativistic quark model 

For sake of simplicity, we present the results for the rates as well the subsequent calcula- 
tions in the rest frame of the decaying particle, although we have also done the calculations 
for a moving one. We sum up our conclusions from the model calculations in the formula 



valid for sufficiently regular potentials. 

A is a typical excitation energy for the mesons levels ; Sing = m h —m c . r indusive and 
its expansion have been calculated as the sum of the widths of any possible state. One can 
see that the terms that are added to 1, correcting the free-quark rate, are corresponding 
to the expected terms of the OPE or more precisely, of the short time expansion which 
plays the role of OPE in this type of models 3 ; this corresponds to the general statement 

3 The operators in the series are local in time, but, unlike in QCD, and not surprisingly for a potential 
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of duality, because duality means basically that the inclusive rate is given by an OPE 
expansion, beginning with the free-quark rate. Therefore, we find duality to be satisfied 
explicitly in the models. Note that there are two small parameters in the OPE expansion, 
one in addition to AQCD/ m b '■ we can choose as the second one ■ 

To present a simpler expansion, we have used above the Shifman-Voloshin (S-V) 
limit [6], where, in QCD, Aqcd <C 8 m Q = m b — m c <^ m b . The second small parameter 
A ££° is then large with respect to kQ CD /m h . In the non relativistic model, the role of 
Aqcd with a light current quark mass is traded for A, the excitation energy with respect 
to the ground state, and the light constituent quark mass m<i- 

The expansion begins with There is no — term. This corresponds in the OPE 

to the fact that we do not find operators with a ^ coefficient. Moreover, there are no 
terms with a positive power n of drriQ, i.e. terms like <y&n ^ > are not present. Explicit 
calculations of the various bound state contributions to T inc i usive in this model show that 
terms of order 0( ^ Sm $ ) (they correspond to form factor effects) and the smaller 0{^r) 

b b 

appear in the exclusive rates, but they cancel in the sum. The respective cancellations 
are ensured respectively by the Bjorken sum rule for 0( ^ Sm $ ) terms and, as we have 

m b 

noticed, by the Voloshin sum rule for 0{^§-) terms [8]. It is found that these sum rules 
are automatically satisfied in the model, which explains the cancellation of the two type 
of terms. 

In the usual 1/rriQ expansion where m c /mb is held fixed, one has instead Sttiq oc itiq, 
and 0{ ^ 5n ffl ) and O(^r-) terms will be counted respectively as (I/itlq) and (1/toq). 

b b 

Since the first type of terms correspond to form factor effects, i.e. bound state effects, they 
are naturally not suppressed by the heavy quark mass in individual state contributions, 
but they must disappear in the total rate, so as to get the free quark rate at leading 
order. On the other hand, the fact that terms of order 0(^g-), which would correspond 

b 

to order 1/rriQ are absent in the ratio R s i, may be expressed as "absence of I/tuq duality 
violation" , because it is in agreement with OPE. As we see, the surviving terms, written 
in eq. (1), are suppressed with respect to the above 0(^t-), by powers of l/m b or by 

b 

powers of A/5rriQ, A being the level spacing, of the order of Aqcd , i-e. they are at most 
of order (l/mg) 2 . 



3 A completely explicit calculation 

To confirm the above statements, it is useful to exhibit the phenomena of concern in 
a completely explicit model. In [11], we have analysed very explicitly the cancellations 
which occur with an harmonic oscillator (H.O.) potential to produce these results. The 
interest of this model is that the truncation of states to the first band of orbital excita- 
tions (lowest D**) becomes exact to the relevant order 1/mjj included, which allows to 
present a very simple analysis ; also it permits to perform completely explicit ana- 
lytical calculations and expansions to high orders or numerical calculations to 
any accuracy. In addition, it is close to the models which are invoked to support the 

theory, their are not local in space. Note that we retain in the paper the term "OPE" , for shortness, but 
also to underline the parallelism with QCD. 
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arguments of [5] , so that one can check with precision the various statements made in this 
paper. In the following is the light quark mass, R 2 , the square of harmonic oscillator 
radius. We find, with A = m ^ R2 in the model (this is the difference between the ground 
state and the first excited level), and with a constant, momentum independent leptonic 
interaction 4 

R s i — 1 H ^— ~ \ ) + smaller terms (2) 

R 2 mf ) \A 5mQ J 

As we said, this expansion can be constructed through an OPE; for example, the first term 
in the bracket corresponds to the kinetic energy operator. In fact we have demonstrated 
that for regular potentials as the H.O. model, the whole OPE series equates the exact 
total width up to terms smaller than any negative power [10]. The conclusion to be 
drawn from eq.(2) is once more that l/mq duality violation, which would show up as a 
larger (9(^r) term, is not present. It remains then to explain where opposite conclusions 

b 

invoking the work [5] may fail. We try to give an answer in sections 5,6. 



4 The condition that many channels be open (SrriQ » 
A) 

But before proposing this answer, we think useful devote a section to an important ques- 
tion about the conditions of validity of the above expansions. Since doubts have indeed 
arisen from apparently contradictory statements, it is useful to discuss in some detail 
why we say that the ratio A/5mq must be small. The smallness of this parameter will 
play a crucial role in the following sections . This condition allows many states to be 
kinematically allowed (which does not imply that transition amplitudes to all of them 
are large) ; this seems to be a natural condition for the OPE expansion, which requires 
inclusiveness, that is a large summation of states. However, the requirement of having 
many channels opened may seem to be contradicted by some well-known results; one may 
think for instance to the well known Shifman and Voloshin result that, in the SV limit, 
the free parton rate equals the sum of rates to ground state mesons D^> only. We show 
that the contradiction is only apparent. 

a) First, let us stress that the SV theorem does not imply that inelastic channels 
are absent ; what happens when many channels are open is that all the contributions, 
other than the one of the ground state, are small with respect to the latter : they include 
a factor (^^-) 2 which is by definition small in the SV limit. 

b) Second, one notes that the Shifman- Voloshin equality of the free quark rate 
and the elastic contribution is violated by a term of same order i^^-) 2 -coming from 

the limited expansion of the form factors. Although (^^-) 2 is small in the SV limit, 
nevertheless, it is a term which is excluded in the totally inclusive rate. And this happens 
because it cancels precisely with the corresponding term from the inelastic contribution, 
leaving only terms small with which are indeed allowed by eq. 1. 

m b 

4 In the paper [11], we have given the results for a more realistic leptonic interaction. The present 
interaction is chosen to fit as much as possible the two-body model of [5] , with a remaining difference as 
regards parameters, which is explained at the end of section 3 
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c) Third, why and when is the the condition bvriQ » A necessary ? The answer 

is : it is in order that the first term in the expansion eq. 1, of the form -\ dominates 

m b 

effectively over the next ones. If 5rriQ is not large with respect to A, one has to take into 
account a whole series of term with increasing powers {j^) n - 

More generally, there are admittedly several numerical situations where the free 
quark rate is saturated very well "by a few channels ", although the mass is not heavy or 
although not many channels are opened. In these latter cases, there is neither an actual 
contradiction with our statement ; what happens is that the difference, although small, 
is no longer indicated by the first term in the above expansion, -^ ; one cannot neglect 

the next terms in the A/5rriQ expansion ; the difference may be practically very small, 
but it comes out as the result of a cancellation of the various orders ; it could be larger, 
or smaller, or of opposite sign with respect to the first term . There is an example below, 
eq.(9), in which we calculate the difference between the exact inclusive rate and the free 
quark decay rate in the realistic situation, where m b — m c is not really small with respect 
to the two masses: then, one observes that the difference is small, but the first term, and 
not even the first two terms in the above S-V series are not sufficient to account for it. 
Other examples, in QCD 2 , given by Lebed and Uraltsev [3], exhibit similar phenomena. 

In short, it is when one desires to use the systematic expansions 1 or 2, that A/Stuq 
must be small, i.e. that the number of open channels must be large. And this expansion 
is the only clearcut, well defined, statement. 



5 Real order of magnitude of the alleged 1/ttiq duality 
violating terms 

Now, let us return to the initial question : what is the explanation of the apparent 
disagreement with the author of [5], who seems to have obtained l/m duality violating 
terms, in contradiction with our SV expansion, where such terms would show up as 
terms and are not found in our calculation ? First, there seems to be a misunderstanding as 
to the magnitude of the effects considered in [5] , induced by the expression "1/ itlq duality 
violation" repeatedly used, and which is very misleading (this section) : the terms are not 
of order but smaller. Second, it seems that in addition very strong cancellations are 
operating in the total rate, which have been missed and lead to a complete disappearance 
of the terms (section 6). We discuss these two points successively. We then proceed to an 
exact calculation which should hopefully discard any remaining doubt (section 7). 

Let us explain our first point. In fact, in [5], it seems not to be actually disputed 
that the OPE expansion is basically right, and that the duality with free quark decay is 
satisfied within the expected accuracy in the region of phase space where the energy release 
is large {t x J^ ax — t 1 / 2 )/A >> 1 (this condition means that many states are kinematically 
allowed). This is certainly true in the physical situation when t is small (i.e. large available 
three-momentum transfer). What may cause problem, according to [5], is only the region 
near t max where this condition is not satisfied. The main point of [5] is that this could 
generate large effects of order 1/mg. However, the expression "order 1/mg" is confusing. 
According to us, the effects which are invoked are of relative order 1 / ijiq only with respect 
to the decay over a small region of phase space : the one where only the ground 
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state can be produced. But -this is the first part of our objection-, these effects must be 
properly referred to the total decay rate, since, in the present debate, it is the totally 
inclusive rate which is to be determined, (the aim is to determine V co by comparison with 
the totally inclusive experimental rate); the debate is about duality violation in the total 
rate ; then we find that the invoked effects are not of order 1/mg with respect to this 
totally inclusive free quark decay rate, but much smaller, by powers of 2A/5mQ, which 
amounts in the standard 1/mg expansion at fixed ratio of heavy masses, to further powers 
of the heavy mass because then, Stuq oc my,. In fact, as we have emphasized above, A/Stuq 
must be taken as small in the general OPE expansion. Also, numerically, the terms are 
small realistically since 2A/<5mg is small (around 2 10 _1 , see the parameters below). This 
we show now. 

Effect I. The first example given by Isgur is that the decrease of the ground state 
contribution with decreasing t (i.e. increasing \q\), due to the form factor, must be 
compensated by the increase of the excited states to maintain duality with free quarks. 
This is exactly guaranteed by the Bjorken sum rule in the heavy quark limit, but it is no 
more exact for finite mass, because there is a region below the D** threshold where only 
the ground state D + D* contribute. Let us write quantitatively the term pointed in [5] , 
with a constant leptonic interaction-the choice of this interaction is not crucial: 

5Tj ~ -K-^ / dt\q\\q\ 2 (3) 



where K is an irrelevant constant, \q\ 2 ~ (5mg) 2 — t ; — p 2 ^- describes the falloff of the 
ground state, and we have approximated the integration limits to the desired accuracy. At 
the lower limit of the integral t = (5mQ — A) 2 , this falloff attains, —p 2 2A ^z Q ', the integral 

is found to be of order with respect to the partial rate over the same region, which 

we can evaluate at leading order by the corresponding free quark rate J^mQ-A) 2 dt\q\ : 

/i- ^ - 3 -^(2A5m Q ) = -l^l (4 ) 

r partial 5 Wlj, 5 TTij, 

where we have used that p 2 A = ^ in the H.O. model. The last expression is seen to be 
of order 1/mg. Indeed, in the 1/mg expansion, <5mg must be counted as (mg) 1 . 

However, this is not what is to be done to evaluate the degree of duality violation 
on the total decay rate ; as we have remarked in our paper [11], the real violation of 
duality is in fact much smaller, because the term must not be divided by the decay rate 
over the above small region of phase space, but rather by the total decay rate, i.e. the 
one over a much larger region of phase space. We end with : 

5T! _ 3 p2 2A5m Q / 2 A \ 3/2 _ 3 m d 5m Q ( 2 A \ 3/2 
F tot 5 m\ \5m Q ) 5 m 2 b \Sm Q J 

Parametrically, this is suppressed with respect to 1/mg, because of the factor (^^) 3 ^ 2 
(5m Q » A). 

Effect II. In another example, relying on a model of two-body decay, the article [5] 
tries to take into account also the effect due to the ma ™ Q terms present in partial rates 

m b 

; such terms, which corresponds to 1/mg, are present separately in the various exclusive 



6 



channels, for instance one has for the ratio Ji^ round state ^ Q f the ground state to the free 
quark decay rates : 

^(ground state) -, 

3m d Sm Q 



b 



but they cancel in the total decay rate. Then the argument would be that if the kine- 
matical situation is such that only the ground state is produced, the total ratio R s i would 
depart from 1 by the term ■2 md6 ™ Q m However, we notice that the effect is for t above the 

b 

D** threshold (mg — m^.*) 2 , i.e. in a limited region of phase space. When we duely refer 
this effect to the total rate, we end with : 



5T H 3 m d 6rriQ f{smQ-A)2 dt\q\ 3 m d 5m Q ( 2A \ 3/2 



\ot 2 m 2 J^ )2 dt\q\ 2 m 2 \5m c 



(7) 



which is once more parametrically smaller than I/tuq, by the factor (J^-) 3 ^ 2 - 

Finally, we conclude that both effects are not C(l/mg), but much smaller. In 
a expansion with fixed ratio of heavy masses, the estimate of these terms, often 

referred to as a potential 1/toq duality violation should be actually much smaller : 



= 0{{l/m Q f 2 ) (8) 



tot 



Numerically too, sticking to the strict harmonic oscillator, and with realistic parameters, 
rrib = 5 GeV, m c = 2 GeV, m d = 0.3 GeV (constituent masses), R 2 = 10 GeV~ 2 , 
giving A = 0.33 GeV (in agreement with heavy light spectroscopy), the effects are small 
: we find respectively : —2.2 10~ 3 and 5.6 10~ 3 . Note that these parameters give a 
low p 2 , much lower than the one observed (p 2 ~ 0.5 instead of p 2 ~ 1) a well known 
defect of the non relativistic approach ; it is then tempting, as done by some authors, 
to introduce additional "fudge" factors to account for a more realistic p 2 , which implies 
that the transition amplitudes will not be given by the basic Hamiltonian ; however such 
a procedure cannot be consistent and will lead to violation of the basic SV sum rules ; 
therefore, although presenting possibly some phenomenological advantage, it is misleading 
for the discussion of duality (in fact, the phenomenological discrepancy for the value of 
p 2 is corrected consistently in a relevant relativistic approach [12], [13]). Note also that 
we have taken the maximum effect by cutting the integrals just at the D** threshold. 



6 Complete cancellation of alleged terms 

Having established the explicit expression of the controversial effects, and demonstrated 
that they are much smaller than stated, one can further wonder whether such terms are 
really present at all in the actual departure from free-quark duality. Indeed, they have 
been picked out as isolated contributions to the decay rate, not as the full contribution 
at some definite order (I/ttiq) 5 ^ 2 . The answer is then that these terms simply do not 
survive when one duely sums up all the contributions of the same order to the 
total rate . Indeed, they correspond to non analytic, fractional powers of (Stuq), and 
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such terms simply do not appear in the expansion, described above, eq. (2), of the ratio 
R s i of the harmonic oscillator model, which derives from a complete and very explicit 
calculation. Moreover, the controversial terms are larger than the last term calculated 
in eq. (2) : they are (l/mf,) 2 l/((5mQ) 1//2 , while we have calculated to (l/mb) 2 l/ (Sttiq) 
included. Therefore, they really cannot be there, they will be cancelled with similar terms 
when duely making the complete calculation of the total rate, as we have done. That such 
a cancellation can occur is already suggested by our calculations of section 5, which shows 
that the two contributions invoked by ref. [5] are of the same type oc !Jh ^¥ L (s^) 3/2 , but 
of opposite sign ; of course, the cancellation is not complete at this stage, but it simply 
means that, in the approach of [5], one has still overlooked some similar effects, which is 
difficult to avoid since one just points out some possible contributions ; one does not rely 
on a systematic enumeration, as we try to do in our expansion. 



7 OPE versus exact numerical calculation at finite 

masses 

Finally, to discard definitively the worry about a possible large (i.e. several 10 -2 ) de- 
parture from OPE, one can ask what is the departure of OPE from the exact rate for 
physical masses and realistic parameters in the model, i.e. far from the S-V limit. An 
exact numerical calculation can be made in the case of the H.O. model. We find that, for 
this H.O. model, close to the one used in [5], but without fudge factors or adjustment of 
phase space, and in an exact calculation where we sum over all the eight allowed H.O. 
levels, and with realistic parameters, R s[ is exceedingly close to the OPE prediction. 
With the above parameters, the departure from 1 is: 

R sl - 1 = -1.33 1(T 3 . (9) 

This already small number has the same order of magnitude as the one predicted by the 
first two terms in the S-V expansion, eq. (2), but the sign is opposite, which suggests the 
importance of the further terms in the S-V expansion ; this is not surprising since mb — m c 
is not sufficiently small to make this expansion very effective : there is no reason in that 
case to reach the correct result with a small number of terms, and one would require a 
calculation of many terms. 

Therefore, for simplicity, instead of the S-V expansion, we now pass to the full 
1/tuq expansion, which we have indeed also performed. Let us recall that in this new 
OPE expansion, it is not assumed that mj — m c is small with respect to m&; each term of 
this new expansion includes a series of many S-V terms. As we have said, we have no term 
of order I/tuq. Considering the l/m^ term in this new expansion, the above discrepancy 
of sign disappears : the OPE estimate of R s i — 1 is now negative and small. The correct 
magnitude itself is reproduced by including also the 1/ttiq term (to the accuracy 10~ 5 ). 
The l/rriq term is very small. Let us insist that in the model, the coefficients of this 
1/rriQ series have been demonstrated to be the ones predicted by the OPE (in fact, the 
short time expansion) , so it means that the possible discrepancy with the OPE series is 
of magnitude 1CT 5 (probably at worst) for realistic masses. 

The fact that the l/mi is very small, while the l/m^ term is of the same order 
as the 1/mg one, is exactly as expected from the arguments of the paper [14] in QCD 
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: in fact, it is the 1/mg term which is abnormally suppressed, whence the abnormally 
small decrease from 1/mg to 1/ttt.q; after that, the terms are decreasing quite rapidly, as 
expected from the effect of inverse powers of mj — m c . 

In summary, in this well-defined and completely calculable quark model, we retrieve 
in an explicit and exact calculation, without any assumption, the conclusion that have 
been inferred in QCD on the OPE of inclusive semi-leptonic decays, and which had been 
controversed on the basis of certain misinterpreted previous quark model calculations . 
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